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What are catamorphisms? 

• A family of functions: often called fold operators. 

• They replace data constructors with arbitrary functions. 

Catamorphisms are useful: 

• they are elegant means of expressing algorithms; 

• they facilitate reasoning about programs; 

• they support program optmization (e.g., deforestation). 
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Goal of this paper 



catamorphisms 
over first-order 
datatypes 




catamorphisms over 
any algebraic type 



Earlier work: Paterson [draft '94], and Meijer & Hutton [FPCA'95] 
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The list catamorphism 

Similar to ML's fold and Haskell's foldr. 
datatype List(cc) = Nil | Cons of a x List(a) 

fun cata(n,c) Nil = n 

I cata(n,c) (Cons(a,r)) = c( a, cata(n,c) r ) 

Function c replaces Cons, value n replaces Nil: 
cata(n,c) (Cons(1,Cons(2,Cons(3 ; Nil)))) = c(1 ,c(2,c(3,n))) 
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Example 

fun filter(p) x = 

cata( Nil, X{a,r). if p(a) then Cons(a,r) else r ) x 

^ N v ' 

n c 



filter(even) (Cons(1 ,Cons(2 ; Cons(3 ; Cons(4,Nil))))) 
= c(1,c(2,c(3,c(4,n)))) 
= Cons(2,Cons(4,Nil)) 
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Catamorphisms can be easily defined for most first-order 
algebraic datatype. 

datatype Term = Const of int 
I Var of string 
I Appl of Term x Term 
I Abs of string x Term 

(Xx. x) 5 ► Appl(Abs("x",Var "x"),Const 5) 



fun cataT(c ; v,p,a) (Const n) = c n 
I cataT(c ; v,p,a) (Var s) = v s 
I cataT(c ; v,p,a) (Appl(x,y)) = p( cataT(c,v,p,a) x, 

cataT^cv^a) y ) 
I cataT(c,v,p,a) (Abs(s,e)) = a( s, cataT(c,v,p,a) e ) 



v 
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Hi gher-order datatypes 

No problem, unless they contain a contravariant recursion: 
datatype T = C of int I D of T -> T 

Why bother? There are 3 examples: 

• higher-order abstract syntax; 

• circular lists; 

• graphs. 
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Hi gher-order abstract syntax 

Pfenning & Elliot [PLDI'88], Wand [POPL'93], Despeyroux & Hirschowitz [LPAR'94] 



Closed terms of ?i-calculus: 

datatype Term = Const of int 

I Appl of Term x Term 
I Abs of Term Term 



(k. x x) 5 ► Appl(Abs(fn x => Appl(x,x)), Const 5) 



v 
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Evaluator over higher-order syntax 

datatype Term = Const of int 

I Appl of Term x Term 
I Abs of Term Term 

datatype Value = Num of int I Fun of Value^Value 

eval: Term Value 

fun eval(Const n) = Num n 

I eval(Appl(f,e))= let val Fun(g) = eval(f) 

in g(eval(e)) end 
I eval(Abs f) = Fun( ... ) 

eval( Appl(Abs(fn x => x), Const 1 ) ) = Num 1 
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Circular Lists 



datatype Clist(oc) = Nil 

I Cons of ax Clist(a) 

I Rec of Clist(a) Clist(a) 

fun head(Cons(a,r)) = a 

| head(Recf) = head(f(Rec f)) 











2 





The circular list [1 ,2,1 ,2,1 ,2,...] is: 

Rec( fn x => Cons(1 ,Cons(2,x))) 

[1,2,3,4,...] is: 

Rec( f n x => Cons( 1 , map(fn z => z+1 ) x ) ) 
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Graphs 

datatype Graph(a) = Node of ax List(Graph(a)) 

| Rec of Graph(a) -» Graph(a) 
| Share of (Graph(a) -> Graph(a)) x Graph(a) 



T(x) 

Share( f n x => T(x), e ) 
e = let val x = e 

in T(x) end x 

0 




Rec( fn x => T(x) ) 
= Y(T) 



fun root(Node(a,r)) = a 
| root(Recf) = root(f(Rec f)) 
| root(Share(f,e)) = root(f e) 

Rec(fn x => Share(fn z => Node(0,[z,Rec(fn y => Node(1 ,[y,z]))]). 

Node(2,[x]))) 



Fegaras & Sheard OGI P0PL96 -11- 



An example of our technique 

• the recursive Term evaluator in the style of Paterson- 
Meijer-Hutton; 

• the problems of this method; 

• the recursive Term evaluator in our style; 

• the catamorphism over Term. 
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The X term evaluator (a la Paterson-Meijer-Hutton) 

datatype Term = Const of int 

I Appl of Term x Term 
I Abs of Term Term 

datatype Value = Num of int I Fun of Value^Value 



fun eval(Const n) = Num n 

I eval(Appl(f,e))= let val Fun(g) = eval(f) 

in g(eval(e)) end 
I eval(Abs f) = Fun( G ) 

Need to define G : Value Value 
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eval(Abs f) = Fun( G ) 



Term 



reify 



Value 



1 ^ Term 
eval 



G 



* Value 



fun eval(Const n) = Num n 

I eval(Appl(f,e))= let val Fun(g) = eval(f) in g(eval(e)) end 
I eval(Abs f) = Fun( eval ° f ° reify ) 

and reify(Num n) = Const n 

I reify(Fun g) = Abs( reify ° g ° eval ) 



v 
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Evaluating (Xx. x) 1 using the eval-reify approach 

eval( Appl(Abs(fn x => x),Const 1) ) 

= let val Fun(g) = eval(Abs(fn x => x)) 
in g(eval(Const 1)) end 

= let val Fun(g) = Fun(fn x => eval(reify(x))) 
in g(eval(Const 1)) end 

= eval(reify( eval(Const 1) )) 

= eval( reify(Num 1) ) 

= eval(Const 1) 

= Num 1 

V ) 



• Double the effort: eval needs reify, print needs parse; 

• Redundancy of computation: 
eval(Appl(Abs f,e)) = eval(f(reify(eval(e)))) 

• reify must be the right inverse of eval: 
Vxg rar?ge(eval): eval(reify(x)) = x 
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Problems of previous work 




Value 



eval(Abs(fn x => x)) 
= Fun( eval ° (f n x => x) ° reify ) 
= Fun( f n x => eval(reify x) ) 
= Fun(fn x => x) 
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Our approach: Using an extra constructor Place to 

shortcut reify 

datatype Term = Const of int 

I Appl of Term x Term 

I Abs of Term Term 

I Place of Value 

fun eval(Const n) = Num n 

I eval(Appl(f,e)) = let val Fun(g) = eval(f) in g(eval(e)) end 
I eval(Abs f) = Fun( eval ° f ° Place ) 
I eval(Place x) = x 

Place satisfies the desired property: eval(Place x) = x 
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Evaluating (Xx. x) 1 using the Place constructor 

eval( Appl(Abs(fn x => x),Const 1) ) 

= let val Fun(g) = eval(Abs(fn x => x)) 
in g(eval(Const 1)) end 

= let val Fun(g) = Fun(fn x => eval(Place(x))) 
in g(eval(Const 1)) end 

= eval(Place( eval(Const 1) )) 

= eval(Place(Num 1)) 

= Num 1 



v 
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The catamorphism over Ter m 

datatype Term(a) = Const of int 

1 Appl of Term(a) x Term(a) 
I Abs of Term(a) -> Term(a) 
I Place of a 

fun cataT(c,p,a) (Const n) = c n 

I cataT(c,p,a) (Appl(x,y))= p( cataT(c,p,a) x, cataT(c,p,a) y ) 
I cataT(c,p,a) (Abs f) = a( cataT(c,p,a) ° f ° Place ) 
I cataT(c,p,a) (Place x) = x 

fun eval x = cataT( Num, fn (Fun(f),e) => f e, Fun ) x 
eval: Term(Value) Value 
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Restrictions 

To avoid problems: 

• The Place constructor should only be used by a cata. 

Solution: make Place and cata primitives and hide 
Place from users. 

• A contravariant variable should not be analyzed. 

eval(Abs(fn x => case x of . . .)) Fun(fn x => eval(case (Place x) of . . .)) 
eval(Abs(fn x => cataT(. . .) x)) -h> Fun(fn x => eval(cataT(. . .)(Place x)) 

Solution: use the type system to detect these cases. 

y 
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Doing this in general: the catamorphism over any 

al gebraic type T 



It is defined in terms of the functor of T. 
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Some background 

The functor for x(oc) maps a function f into the function |[x(a)](f): 
f: t!-»t2 ^ Ix(a)](f): Tft)-^ 



M(f) 


= f 


Notation: 


lT 1 XT 2 ](f) 


= M(f) X [T 2 ](f) 


fxg = A<x,y). ( f x, g y ) 


U(T)](f) 


= map T ([T](f)) 


id = Xx. x 


M(f) 


= id 





e.g., |[ int x a x list(a) ](f) = id x f x map(f) 

= A,(x,y,z). ( x, f(y), map(f) z ) 



Properties: |x]](id) = id 



IxI(f) 0 IxI(h) = [x](f°h) 
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Extending functors to arrow types 



Need to differentiate positive from negative types: 
+ - + 
(int — > bool) — > int 



Positive and negative instances of a type variable a: 

+ - - + 
type T(a) = (a^a)^(a^a) 

X v ' « v ' 

+ 
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Bifunctors 

The bifunctor of a type t(cc): 

f: t!-»t2 A g: t 2 ->tj => [T(a)](f,g): x(ti)^T(t 2 ) 



M(f,g) 


f 


lT 1 XT 2 ](f,g) = 


[Tll(f,g) X [T 2 l(f,g) 


IXl^T 2 l(f.g) = 


Xh. [T 2 l(f,g) ° h ° [Ti](g,f) 


[T(T)](f,g) = 


map T ([x](f,g)) 


M(f,g) 


id 


Properties: 




M(id,id) 


= id 


M(f,g) ° Ix](h,k) 


= [x](f 0 h,k 0 g) 
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The catamorphism over a type T 

Type T has n data constructors T|(T)-»T, 
plus one extra constructor Place : 



cata T (f 1? . 


f ) 0 C = 


fi ° [Ti(a)]( cata T (f 1 ,...,f n ), Place 1 ) 


cata T (f l5 . 


.,f n )° Place 1 = 


id 



datatype T(P) = C of ( T(P) ^ T(P) ) ^ T(P) 
I PlaceTof(3 

fun cataT(f)(C g) = f(cataT(f) ° g ° (fn h => PlaceT ° h ° cataT(f))) 
| cataT(f) (PlaceT x) = x 
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Conclusion 

Simulating the right inverse of a function with an extra 
data constructor Place is useful for: 

• meta-programming; 

• circular lists and graphs. 

It might also be useful for: 

• higher-order unification; 

• fusing programs using parametricity laws. 

' 
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Equality between terms 



fun eqT(Const(n),Const(m)) k = ( n=m ) 

| eqT(Appl(a,b) I Appl(a\b') k = eqT(a,a') k A eqT(b,b') k 

| eqT(Abs(f),Abs(g)) k = eqT(f(Place k) ; g(Place k)) (k+1 ) 

| eqT(Place n, Place m) k = ( n=m ) 

| eqT__ = false 



eqT( Abs(fn x => Appl(x,x)), Abs(fn y => Appl(y,y)) ) 0 

= eqT( Appl(Place 0, Place 0), Appl(Place 0, Place 0) ) 1 
= eqT(Place 0, Place 0) 1 A eqT(Place 0, Place 0) 1 
= (0=0) A (0=0) 
= true 
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Program fusion 

The Y combinator for functions: 
Y(f) x = f(Y(f)) x 

Y: Voc,p\ ((a^(3)^a^(3)^a^(3 

The parametricity theorem for Y is the unfold- simplify-f old law: 
Vg,U,h: g°(f h) = ^ (g 0 h) ^ g°(Yf) = Y^ 

To find a solution for (j): 

1) choose Place such that g ° Place = id, 

2) choose k such that h = Place ° k, then: 

g ° (Y f) = Y(^k. g ° (f (Place ° k))) 
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datatype list(a,p) = Nil | Cons of a x list(a,(3) | Place of (3 

fun lenNil =0 fun inc Nil = Nil 

I len (Cons(a,r)) =l+(lenx) I inc (Cons(a,r)) = Cons(a+l,inc r) 

I len (Place x) = x I inc (Place x) = x 



f x = len(inc x) 

fun f Nil =len(Nil) 

I f (Cons(a,r)) = len( Cons(a+ 1 ,Place(f r)) ) 

I f (Place x) = x 

f [ai,...,a,J creates n garbage Cons-Place cells: 

Cons( a i? Place(n-i) ) 



Fegaras & Sheard OGI POPL96 _ 29 



Solving parametricity equations 

The parametricity theorem for the type Vcc. x(a) that 
relates x and y is: 

x = [x(a)]( f, Place ) y 
f ° Place = id 

Example: 

g: Voc: x 1 ^x 2 ^...^x n ^a 

f(g Xl ... x n ) = g ([XiK f, Place ) Xl ) ... ([xj( f, Place ) x n ) 

cata: ax((3xa^a)^list((3)^oc 

f(cata(n,c) x) = cata( f(n), ?i(a,r). f(c(a,Place r)) ) x 



Fegaras & Sheard OGI POPL96 . 3Q 



